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ABSTRACT: We investigate a detector scheme designed to measure the arrival of 
a particle at x = during a finite time interval. The detector consists of a two state 
system which undergoes a transition from one state to the other when the particle crosses 
X = 0, and possesses the realistic feature that it is effectively irreversible as a result of 
being coupled to a large environment. The probabilities for crossing or not crossing x = 
thereby derived coincide with earlier phenomenologically proposed expressions involving 
a complex potential. The probabilities are compared with similar previously proposed 
expressions involving sums over paths, and a connection with time operator approaches is 
also indicated. 



1. INTRODUCTION 



An enduring class of questions in non-relativistic quantum theory are those that involve 
time in a non-trivial way [1,2,3,4]. Of these, the question of tunneling time is perhaps the 
most important [5,6]. But some of the basic issues are most simply seen through the 
question, what is the probability that a particle enters a region of space for the first time 
during a given time interval? What makes this sort of problem interesting is that standard 
quantum mechanics seems to handle it with some difficulty and there does not appear to 
be a unique answer - classically equivalent methods of assigning the arrival time can differ 
at the quantum level. 

One of the key sources of difficulty with defining arrival times in quantum theory is that 
they involve specification of positions at two moments of time: to state that the particle 
enters a given spatial region at time t means that it is outside the region immediately 
before t and inside it immediately after t. Since positions at different moments of time do 
not commute, we do not expect to be able to associate a single hermitian operator with 
the arrival time. Of course, many different methods of defining arrival time then naturally 
suggest themselves, but, like the question of specifying phase space locations in quantum 
mechanics, they are not all equivalent. 

To be specific, in this paper we consider the following question: Suppose we have an 
initial state with support only in a; > 0. Then what is the probability that the particle 
enters the region a; < at any time during the time inter al [0, r]? 

A number of previous papers have addressed this particular problem using path inte- 
grals [7,8,9]. (These approaches are closely related to the decoherent histories approach to 
quantum theory [10,11]). The amplitudes for crossing and not crossing x = are obtained 
by summing over paths which either always enter or never enter a; < 0, and probabilities 
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are then obtained by squaring amplitudes in the normal way. However, due to interference 
between paths, the resultant probabilities do not add up to 1, so cannot be regarded as 
true probabilities. A way past this difficulty was explored in Ref.[12]. There, the point 
particle system was coupled to a thermal environment to induce decoherence of different 
paths in configuration space, and the correct probability sum rules were restored. Although 
this approach produced mathematically viable candidates for the probabilities of crossing 
and not crossing, they depend to some degree on the environment, and it is by no means 
clear how the results correspond to a particular type of measurement, even an idealized 
one. General theorems exist showing that decoherence of histories implies the existence 
at the final end of the histories of a record storing the information about the decohered 
histories [10], but these have been explicitly found only in a few simple cases (see Ref.[13], 
for example). For these reasons, it is of interest to compare these path integral approaches 
with a completely different approach involving a specific model of a detector. 

Let us therefore introduce a model detector which is coupled to the particle in the 
region x < 0, and such that it undergoes a transition when the coupling is switched on. 
Such an approach has certainly been considered before (see, e.g., Ref.[4]). The particle 
could, for example, be coupled to a simple two-level system that flips from one level to 
the other when the particle is detected. One of the difficulites of many detector models, 
however, is that if they are modeled by unitary quantum mechanics, the possibility of 
the reverse transition exists. Because quantum mechanics is fundamentally reversible, 
the detector could return to the undetected state under its self-dynamics, even when the 
particle has interacted with it. 

To get around this difficulty, we appeal to the fact that realistic detectors have a very 
large number of degrees of freedom, and are therefore effectively irreversible. They are 
designed so that there is an overwhelming large probability for them to make a transition 
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in one direction rather than its reverse. In this paper we introduce a simple model detector 
that has this property. This is achieved by coupling a two-level system detector to a large 
environment, which makes its evolution effectively irreversible. The description of this 
system is easily obtained using some standard machinery of open quantum systems, and 
the resulting master equation for the particle coupled to the detector actually has the 
Lindblad form [14]. 

We are not concerned with a specific experimental arrangement, but rather, as is com- 
mon in quantum measurement theory, an idealized model which has as many physically 
realistic features one can reasonably incorporate. In particular, in contrast to most mea- 
surement models discussed in the literature, it has the key property of irreversibility. 

The detector is described in Section II. On solving the detector dynamics, an expres- 
sion for the probability of entering a spacetime region is obtained. It has the appearance 
of a probability obtained from a wave function satisfying a Schrodinger equation with an 
imaginary contribution to the potential, which has previously been proposed as a phe- 
nomenological device [1,15,16]. Our detector model therefore justifies previously used 
phenomenological approaches. 

The probabilities obtained are also readily compared with the results of the path 
integral approaches, and the comparison sheds some light on the shortcomings of the latter. 
This is described in Section III. A brief mention is also made of the possible connection 
with time operators. 



2. THE DETECTOR MODEL 



In this section we describe a detector designed to make a permanent record of whether 
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or not a particle enters the region a; < at any time during a given finite time interval. We 
take the detector to be a two- level system, with levels |1) and |0), representing the states 
of no detection and detection, respectively. Introduce the raising and lowering operators 

a+=|l)(0|, a-_ = |0)(l| (2.1) 

and let the Hamiltonian of the detector be = ^Tiujazi where 

a,= |l)(l|-|0)(0| (2.2) 

so |0) and |1) are eigenstates of with eigenvalues —\fujj and ^Tiuo resepectively. We 
would like to couple the detector to a free particle in such a way that the detector makes 
an essentially irreversible transition from |1) to |0) if the particle enters a; < 0, and remains 
in |1) otherwise. This can be arranged by coupling the detector to a large environment 
of oscillators in their ground state, with a coupling proportional to 9{—x). This means 
that if the particle enters the region a; < 0, the detector becomes coupled to the large 
environment causing it to undergo a transition. Since the environment is in its ground 
state, if the detector initial state is the higher energy state |1) it will, with overwhelming 
probability, make a transition from |1) to the lower energy state |0). A possible Hamiltonian 
describing this process for the three-component system is 

H = Hs + Ha + Hs + V{x)Hds (2.3) 

where the first three terms are the Hamiltonians of the particle, detector and environment 
respectively, and is the interaction Hamiltonian of the detector and its environment. 
The simplest choice of environment is a collection of harmonic oscillators, 

Hg = ^ TiuJnalan (2.4) 
n 
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and we take the coupling to the detector to be via the interaction 

Hd£ = X] ^ («n^-ai + KnO+On) (2.5) 

n 

An environment consisting of an electromagnetic field, for example, would give terms of 
this general form. V{x) is a potential concentrated in a; < (and we will eventually make 
the simplest choice, V{x) = d{—x), but for the moment we keep it more general). The 
important feature is that the interaction between the detector and its environment, causing 
the detector to undergo a transition, is switched on only when the particle is in a; < 0. 

A similar although more elaborate model particle detector has been previously studied 
by Schulman [17] (see also Refs.[18]). It consists of a lattice of spins in a metastable state 
interacting with the lattice's phonon modes, with interactions essentially the same as in 
Eq.(2.3). His approach has the advantage that the amplification of a microscopic event 
as well as the irreversibility of the measurement is modeled, but it does not appear to be 
possible to solve it as explicitly as the simpler model considered here. For an even more 
elaborate model see Ref.[19]. 

We are interested in the reduced dynamics of the particle and detector with the en- 
vironment traced out. Hence we seek a master equation for the reduced density operator 
p of the particle and detector. With the above choices for and H^£, the derivation of 
the master equation is standard [20,21] and will not be repeated here. There is the small 
complication of the factor of V{x) in the interaction term, but this is readily accommo- 
dated. We assume a factored initial state, and we assume that the environment starts 
out in the ground state. In a Markovian approximation (essentially the assumption that 
the environment dynamics is much faster than detector or particle dynamics), and in the 
approximation of weak detector-environment coupling, the master equation is 

p=-'-[Hs + Hd,p] - I {y'^{x)a+a-p + pa+a-V'^{x) - 2V{x)a-pa+V{x)) (2.6) 
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Here, 7 is a phenemonological constant determined by the distribution of oscillators in the 
environment and underlying coupling constants. The frequency u) in is also renormal- 
ized to a new value a;'. 

Eq.(2.6) is the sought-after description of a particle coupled to an effectively irreversible 
detector in the region a; < 0. In the dynamics of the detector plus environment only (i.e., 
with V = 1 and Hs = 0), it is readily shown that every initial state tends to the state 
|0)(0| on a timescale 7""*^. With the particle coupled in, as in (2.6), if the initial state of 
the detector is chosen to be it undergoes an irreversible transition to the state |0)(0| 

if the particle enters x < 0, and remains in its initial state otherwise. 

Although we have outlined the derivation of this master equation for a particular 
choice of environment and detector-environment coupling, we expect that the form of the 
equation is more general. It is well-known that, after tracing out the environment, and in 
the approximation that the evolution is Markovian, the reduced density operator p of the 
particle and detector must evolve according to a master equation of the Lindblad form 

P = -^[^s + H^, p] + Y^ (impLln - ^L^nLmP - ^pLjnLmj (2.7) 

This is the most general Markovian evolution equation preserving the positivity, hermiticity 
and trace of p. The operators Lm model the effects of the environment [14]. The form 
(2.6) is recovered with a single Lindblad operator L — 'y2V{x)a-. A similar detection 
scheme based on a postulated master equation similar to (2.7) was previously considered 
in Ref. [22,23], although the resultant expressions for arrival time probability given below 
were not derived (and no microscopic origin of the equation was given). 

Eq.(2.6) is easily solved by writing 

p = PU ® |1)(1| + poi ® |0)(1| + pio ® |1)(0| + POO ® |0)(0| (2.8) 
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where 



Pll = 


* rtj- 1 
-^[-f^s,PllJ 


1 / -r r / \ . x7-/\\ 

-^{V{x)pii + piiV{x)) 


(2.9) 


POl = 


i 


7 , .Tiuj' 

^pQlV{x)+l 2 POl 


(2.10) 


PlO = 


-^[^s,Plo] 


^V{x)piQ I 2 PlO 


(2.11) 


POO = 


i 

-^[-f^s,P00j 


+ lV{x)piiV{x) 


(2.12) 



(where we have now set V{x) = 9{—x), so that V = V). We suppose that the particle 
starts out in an initial state \^o)-, hence the above equations are to be solved subject to 
the initial condition, 

p(0) = |*o)(*ol®|l)(l| (2.13) 
The probability of finding the detector in the unregistered state |1) at time r is 

/oo 
dx pii{x,x,t) (2.14) 
-oo 

and the probability of finding it in the registered state |0), is 

/oo 
dx poo{x,x,t) (2.15) 
-00 

(where the trace is over the particle Hilbert space). Clearly + — 1, since Trp = 1. 

Note that the probability for no detection includes an integral over x < and pii(a^, x, r) 
is not necessarily zero for a; < 0. There is therefore the possibility that the particle could 
enter the region a; < without the detector registering the fact. This is however, to be 
expected of a realistic detector - there is some probability that it will fail to do what it is 
supposed to do. The probability of this happening is typically small, and indeed, compu- 
tation of this probability provides a useful check on the efficiency of the detector (although 
below we will check detector efficiency in a different way). 
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The formal solution to Eq.(2.9) for pn may be written 



pn{t) = exp i^—Hst - |wj pii(O) exp (^^Hst - 'Ivtj (2.16) 

What is particularly interesting about this expression is that it can be factored into a pure 
state. Let pu = |^')(*|. Then, noting that pii(O) = |^'o)(*ol7 Eq.(2.16) is equivalent to 



\^{t))^e^p[~Hst-^Vt]\^o) 



(2.17) 



The probability for no detection is then 



Pnd 



oo 



dx |*(a;,r)p 



(2.18) 



— oo 



The pure state (2.17) evolves according to a Schrodinger equation with an imaginary 
contribution to the potential, —^ih^V. Complex potentials have been used previously in 
this context, as phenomenological devices, to imitate absorbing boundary conditions (see, 
for example Refs. [1,15,16]). Here, the appearance of a complex potential is derived from 
the master equation of a particle coupled to an irreversible detector, which in turn may be 
derived from the unitary dynamics of the combined particle-detector-environment system. 

Eq.(2.12) for poo ni^y be formally solved to yield 

pm{t)=l f^dt'e^p (^~Hs{t-t')^ V{x)pii{t')V{x)e^p (^-Hs{t-t')^ (2.19) 

(recalling that poo(O) = 0)- Inserting the solution for pii{t'), the probability for detection 
may be written, 
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Pd = l 



dx 



Jo J- 

rT rO 

^ dt dx 
Jo J-oo 

^ dt dx\^(x,t)f 
Jo J-oo 



(x|exp ( -^Hs{T-t)j V{x)exp i-'-Hst-^Vt ) |*o) 



[x\e^p[~Hst-^Vt]\^Q) 



(2.20) 



where '^{x,t) is the wave function (2.17). The expression for the probability for detection 
has an appeahng form: it is the integral of |\E'(a;, t) p over the space-time region of interest. 
It is crucially important, however, that the wave function satisfies not the usual Schr5dinger 
equation, but one with an imaginary contribution to the potential modeling the detector. 

It is useful to write the probabilities for detection and no detection in the form, 

Prf = IY(n |*o)(*o|), Pnd = ^(^^ l*0)(*0|) (2.21) 

where 



n = J dt exp y-^Hst - -^Vtj -fV exp y-^Hst - jA^tj (2.22) 
= exp (^^HsT ^Vr^ exp (^-^H,t ^ ^Vr^ (2.23) 

The first of these expressions follows from (2.20) from the properties of the trace, and 
using the fact that V"^ = V. 

Q and Q are clearly not projection operators, although their properties are close to 
those of projectors. They are both positive operators and + = 1. The latter follows 
by integrating the identity 

f = - exp (jrH.r - IVr^ ,V exp [-'-H.r - Ivr) (2.24) 

Moreover, these operators clearly have the desired localization properties on histories of 
particle positions, as will be seen most clearly in the path integral expression of the next 
section. We do not expect to be able to associate a true projection operator with the 
arrival time, but here we have found a POVM, which is the next best thing. 

Eqs.(2.21)-(2.23) are the main result of this section: expressions for the probabilities 
of entering or not entering a region of spacetime, derived from an irreversible detector 
model. 
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We now consider the issue of the efficiency of the detector. A simple way to do this is 
to introduce a second detector identical to the first and located in the region a; > 0. Since 
the entire x-aixis is now monitored, the probability that neither detector registers during 
the time interval is then a measure of the degree to which the detector will fail. 

With two detectors in place, the master equation now is, 

i 'y 

p = --[Hs,p]- -{9{-x)a+a-p + pa+a-9{-x) - 29{-x)a-pa+9{-x)) 

-^{9{x)a+a-p + pa+a-9{x) - 29{x)a-pa+9{x)) (2.25) 

where a+, d- are the raising and lowering operators for the detector in a; > 0. This equation 
is solved like (2.6), by writing 

P = Pnd^|l)(l|^|l)(l| + --- (2.26) 

We are interested only in the probability that neither detector registers, so we omit the 
explicit form of the other terms in (2.26). It is readily shown that 

Pnd = -f\^s, Pnd] - IPnd (2-27) 

The probability of no detection is Trp^^^^, and from (2.27), it clearly decays like e~'^*. Hence 
the detector functions efficiently if the total time duration r is much greater than 7^^. 
(The potential inefficiency of the detector for r not sufficiently large compared to 7"-*^ 
corresponds to the difficulties in defining the time operator at low momenta [3].) 

The evolution according to (2.17) for the case in which V{x) is a real step function was 
studied in detail by AUcock [1] , who was consequently pessimistic about the possibility of 
defining arrival time. This is partly because 7 needs to be large for accurate detection, 
but in this case reflection from the potential is high and not all of the incoming flux 
is absorbed. Subsequent authors have shown that his pessimism was misplaced, if more 
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general potentials V{x) are permitted, which are smoother and may also need to be complex 
(see, for example, Refs. [15,24]). The considerations of this paper readily generalize to this 
case by replacing Eq.(2.6) with the equation 

p = -^[Hs,p] - I (v'fVa+a-p + pa+a-V^^V - 2Va-pa+v'f^ (2.28) 

which is still of the Lindblad form, but V{x) is generally a non-hermitian operator. It is 
less clear what sort of detector-environment coupling this corresponds to in Eq.(2.3) and 
this is certainly of interest to investigate. Note also that ^ V, but the above results 
are readily modified to incorporate this. 

The detector described here measures whether the particle entered the region a; < 
at any time during [0,t], for r >> 7~^. It could easily be extended to give more precise 
information about the time at which the particle enters a; < by using a series of similar 
detectors, but with a time-dependent coupling to the environment, so that the detectors 
can be switched on or off at a succession of times ti, t2, t^ - ■ ■ (separated in time by at least 
7"-*^) in the interval [0,t]. 



3. COMPARISON WITH OTHER APPROACHES 

The expressions we have derived for detection and no detection bear a close resemblance 
to previously derived path integral expressions for the probabilities of entering or not 
entering the region a; < 0, so we now carry out a comparison. 

For simplicity let the initial and final points lie in a; > 0. The amplitude for remaining 
restricted to the region a; > is 

gr{xf,T\xo, 0) = J^Vx{t)exi)(^^S[x{t)]^ (3.1) 
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where the sum is over all paths in a; > 0. The amplitude for crossing a; = at some time 
in the interval [0, r] is 

gc{xf,T\xo,0) = Jvxit)e^p(^^S[x{t)]^ (3.2) 

where the sum is over paths that spend some time in x < (see Refs. [25,26,27] for details 
of construction of these objects). Clearly gr + dc = 9, where g denotes the unrestricted 
propagator, given by a sum over all paths. The probabilities for crossing and not crossing 
X = are then obtained from these propagators, by attaching an initial state, squaring the 
amplitudes, and then summing over final positions in the usual way. However, as stated 
earlier, the resultant probabilities for crossing and not crossing computed in this way do 
not sum to 1, because of interference between the different types of paths. 

Now we compare with the measurement model of the previous section. The probabil- 
ities computed here for detection or no detection in the region x < automatically sum 
to 1. The probability for no detection may computed from (2.17). The evolution operator 
that appears there may be written in path integral form, 

(i T \ 

= J Vx{t) exp (^^S[x{t)] -\f^dt V{x{t))^ (3.3) 

The sum is over all paths x{t) connecting xq at time to at time r. But it is clear 
that the potential V{x) suppresses contributions from paths that enter a; < 0. Split 
the paths summed over into the two classes r and c, as above. (For simplicity, we take 
a;o > 0, a;y > 0). Noting that F = in a; > 0, the path integral becomes. 



S'nd(a^/>Tko,0) = j Vx{t)exp (^^S[x{t)] 
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Comparing with Eq.(3.1), we see that (3.4) differs from it by the presence of the second 
term. In the second term, every path in the sum has a section lying in the region a: < and 
an exponential suppression factor will come into play, g^^ and gr exactly coincide in the 
limit 7^00. The resultant probabilities are, however, not very interesting [7]. Furthermore, 
as stated, large 7 means that most of the incoming wave packets are reflected rather than 
absorbed by the detector. This means that the second term in Eq.(3.4) will generally 
be significant and g^d and gr will not be close. From this we conclude that the purely 
path integral approaches to defining the arrival time, as expressed by (3.1) and (3.2), 
are actually rather removed from the more physically motivated expressions using a model 
detector derived in this paper. Including a physical mechanism for decoherence in the path 
integral approach [12] yields more sensible results, but they are not very closely related to 
the detector results (a comparison is carried out in an earlier version of this paper [27]). 

We may also write the POVM for no detection, (2.23), in another more enlightening 
form. Note that 



Here, T denotes time ordering, and xt = x + pt/m is the position operator at time t in 
the Heisenberg picture. Splitting the time interval into small units 6t, U is therefore a 
time-ordered product of operators of the form ex'p{—jStV/2). But with the choice used 
here, V{x) = d{—x), we have = V so 




(3.5) 



where 




(3.6) 




e{x) + e{-x)e-^^^l'^ 



(3.7) 
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This is therefore a projector onto the positive a;-axis plus an exponentially smaller projector 
onto the negative a;-axis. Naive expectations would lead one to guess the first term, but the 
addition of the second one appears to be important in making the probabilities well-defined 
(see Ref.[26] for related formulae). 

From the above we see that the probabilities for detection and no detection depend on 
X and p only through the operators 0{xt) at a series of times. This leads us to a connection 
with time operators. For classical quantities x, p, we have 

0(a: + ^) = 9{— + t) (3.8) 
m p 

(for x,p > 0). The point is that the quantity mx/p is precisely the classical arrival time. 
It is the quantity that numerous authors have, not without serious difficulties, attempted 
to elevate to the status of a quantum operator to define quantum arrival times (see, for 
example, Ref.[3]). The connection between the expressions derived here and those derived 
using the time operator may therefore be obtained by investigating the extent to which 
the classical relation (3.8) persists in the quantum theory. This will be pursued elsewhere. 



4. SUMMARY AND CONCLUSIONS 

We have presented a detector model for the measurement of arrival time in quantum 
theory. It possesses the realistic feature of being effectively irreversible. The results of 
the scheme connect very nicely with previous approaches involving a postulated complex 
potential to imitate the effects of a detector. 

The detector model was compared with previous approaches involving sums over paths, 
and the detector model exposes the limitations of the latter. Some indication of the possible 
connection with time operators was also given. 
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